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Abstract. Using the Feynman parameter method, we have calculated in an elegant manner a set of 
one— loop box scalar integrals with massless internal lines, but containing 0, 1, 2, or 3 external massive 
lines. To treat IR divergences (both soft and collinear) , the dimensional regularization method has been em- 
ployed. The results for these integrals, which appear in the process of evaluating one— loop (N > 5)— point 
integrals and in subdiagrams in QCD loop calculations, have been obtained for arbitrary values of the 
relevant kinematic variables and presented in a simple and compact form. 



1 Introduction 

Scattering processes are one of the most important sources 
of information on short— distance physics and have played 
a vital role in establishing the fundamental interactions of 
nature. In testing various aspects of QCD, the scattering 
processes in which the total number of particles in the 
initial and final states is N > 5 (like 2 — > 3, 2 — > 4, etc.) 
are becoming increasingly important. 

The techniques for calculating the tree level ampli- 
tudes involving a large number of particles in the final 
state are well established [EJ. Owing to the well— known 
fact that the LO predictions in perturbative QCD do not 
have much predictive power, the inclusion of higher— order 
corrections is essential for many reasons. In general, higher- 
order corrections have a stabilizing effect reducing the de- 
pendence of the LO predictions on the renormalization 
and factorization scales and the renormalization scheme. 
Therefore, to achieve a complete confrontation between 
theoretical predictions and experimental data, it is very 
important to know the size of radiative corrections to the 
LO predictions. 

Obtaining radiative corrections requires evaluation of 
one— loop integrals which arise from a Feynman diagra- 
matic approach. The case of massless internal lines is of 
special interest, because we often deal with either really 
massless particles (gluons) or particles whose masses can 
be neglected in high— energy processes (quarks). The main 
technical difficulty in obtaining the NLO corrections con- 
sists in the treatment of the occurring TV— point tensor and 
scalar integrals with massless internal lines. In QCD, ten- 
sor integrals appear in which the N— point integral may 
contain up to N powers of the loop momentum in the 
numerator of the integrand. Since these integrals contain 
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IR divergences, they need to be calculated in an arbitrary 
number of dimensions and the standard methods of |2| 
cannot be directly applied. 

Various approaches have been proposed for reducing 
the dimensionally regulated (N > 5)— point integrals to 
a linear combination of N— and lower— point scalar inte- 
grals multiplied by tensor structures made from the metric 
tensor g^ v and external momenta f|,|],[||(].0- 

It has also been shown that the general (N > 5)— point 
scalar one— loop integral can be recursively represented as 
a cyclically symmetric combination of (N — 1)— point in- 
tegrals, provided the external momenta are kept in four 
dimensions @,||,[|. Consequently, all scalar integrals oc- 
curring in the computation of an arbitrary one— loop (N > 
5)— point integral with massless internal lines can be re- 
duced to a sum over a set of basic scalar box (N = 4) 
integrals with rational coefficients depending on the ex- 
ternal momenta and the dimensionality of space— time. 
This set of diagrams includes IR divergent box integrals 
with massless internal lines but containing 0, 1, 2, and 3 
external masses and the IR finite box integral with four 
external masses. 

The IR finite box integral has been evaluated in Ref. 
M, and written in a more compact form in Ref. ||. The 
results for the IR divergent box integrals with no external 
masses and with one external mass have been obtained in 
Rcfs. |l(],[llj • All IR divergent box integrals have been con- 
sidered in Ref. [[l2| , using the partial differential equation 
technique. However, the results obtained for the integrals 
containing two and three external masses are strictly cor- 
rect only in the Euclidean region where all relevant kine- 
matical variables are negative. 

Being very nontrivial, impossible to check numerically, 
and of fundamental importance for one— loop calculations 
in perturbative QCD with massless quarks, it is absolutely 
essential that these integrals should be evaluated and the 
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results of Ref. |f^] should be checked using independent 
techniques. 

In this paper we recalculate the IR one— loop box scalar 
integrals in an elegant manner using dimensional regu- 
larisation and the Feynman parameter method, and give 
results in a simple and compact form. A characteristic 
feature of our calculation is that the causal ie has been 
systematically kept throughout the calculation, so that 
the results obtained are valid for arbitrary values of the 
relevant kinematic variables. 

The paper is organized as follows. Section 2 is de- 
voted to introducing the notation and to some prelimi- 
nary considerations. In Sec. 3, using the Feynman param- 
eter method and the dimensional regularization method, 
we evaluate the IR divergent one— loop box Feynman in- 
tegrals with massless internal lines but containing 0, 1, 2, 
and 3 massive external lines, and compare our results with 
the corresponding ones obtained in Ref. 12 . Section 4 is 
devoted to some concluding remarks. An analytical proof 
of the equivalence of our results to those obtained in Ref. 
]l2| for Euclidean kinematics is given in Appendix A. For 
the reader's convenience, in Appendix B we present the 
closed— form expressions for the IR divergent one— loop 
scalar box integrals evaluated in this paper, i.e., with all 
poles in em ; = D/2 — 2 manifest, and with all functions of 
the kinematic variables expressed in terms of logarithms 
and dilogarithms. 



2 Preliminaries 

The massless scalar one— loop box integral in D— dimen- 
sional space— time is given by 

where pi, 7=1,2,3,4 are the external momenta, I is the loop 
momentum, and /i is the usual dimensional regularization 
scale. As indicated in Fig. [j], all external momenta are 
taken to be incoming, so that the massless propagators 
have the form 

Ai = I 2 + ie, 

A 2 = {l+pif + ie, 

A 3 = {I+P1+P2) 2 +ie, 

A 4 = (l+Pi +P2 +Ps) 2 +ic (2) 

Combining the denominators with the help of the Feyn- 
man parametrization formula 



1 



A 1 A 2 A 3 A i 

dxida^da^da^ , , , , 

10 (x 1 A 1 +x 2 A 2 +x 3 A 3 +x 4: A 4 ) 4 ' 



(3) 



3! S(xi + x 2 + x 3 + xa - 1) 



performing the D— dimensional loop momentum integra- 
tion using 

r d D i 



1 



T(4-D/2) 



(2w) D (I 2 - M 2 + ie) 4 (4ir) D / 2 3! (M 2 - ie) 4 -- / 2 ' 

(4) 



Pi 



P'2 



l+Pl 




1+PI+P2+P3 



Pi 



P3 



Fig. 1. Basic one- loop box diagram. 

introducing the external "masses" 

p 2 = m 2 (» = 1,2,3,4), 
and the Mandelstam variables 

s = (Pi +P2) 2 , t=(p 2 +p 3 ) 2 , 



(5) 



(G) 



one readily finds that the scalar integral in ([[]) can be 
written in the form 

T( , 2 2 2 i ^(4 - D/2) 

dxida;2da;3da;4 5{x\ + X2 + x 3 + x 4 — 1) 

x (— xix 3 s — X2X4 t — x\X2 rn\ 

2 2 2 • N-D/2-4 

—x 2 x 3 m 2 — x 3 x 4 m 3 — X\X 4 m 4 — ie] . (7) 

This is the basic four— point "scalar" parametric integral, 
serving as a starting point for our further considerations. 

Depending on the number of the external massless 
lines, we distinguish six special cases of this integral. Fol- 
lowing the notation of Ref. , we denote these integrals 

by 



r4m 
J 4 


j / 1 2 2 2 2\ 

= l4(s,t;m 1 ,m 2 ,m 3 ,m 4 ), 


(8) 


r3m 
1 A 


= l 4 (s,t;0,ml,ml,ml), 


(9) 


j2mh 
M 


= I 4 (s,t; 0,0, 7713,7714), 


(10) 


r2me 
J 4 


= h(s,t; 0,7773, 0,7774), 


(11) 


rim 
1 A 


= h(s,t; 0,0,0, ml), 


(12) 


rOm 
i 4 


= i4(M;0,0,0,0,), 


(13) 



and refer to them as the four— mass scalar box integral, 
the three— mass box integral, the "hard" two— mass box 
integral, the "easy" two— mass box integral, the one— mass 
box integral, and the massless box integral, respectively. 

These six box integrals constitute the fundamental set 
of integrals, in the sense that an arbitrary one— loop N(> 
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5)— point integral with massless internal lines can be repre- 
sented in a unique way as a linear combination of these in- 
tegrals with the coefficients being rational functions of the 
relevant kinematic variables and the number of space— di- 
mensions D. 

These integrals arise from the Feynman diagrams de- 
picted in Fig. H, formally corresponding to the scalar mass- 
less ^ 3 theory. The thick lines in these diagrams denote 
the massive (off— shell) external lines. As it is seen from 
Fig. ||, there are two distinct configurations related to the 
case when two external lines are massless: the adjacent 



= m| = 0) and the opposite box dia- 
0). They correspond to the hard and 



box diagram (m\ 
gram (mf = m\ = 
easy two— mass box integrals, respectively. 

When evaluating the diagrams of Fig. g, one comes 
across IR singularities (both collinear and soft). Let us re- 
call the circumstances under which IR singularities appear 
in a Feynman diagram. When an on— shell quark of mo- 
mentum p emits a gluon of momentum k, then IR singu- 
larities appear as a result of vanishing of the quark prop- 
agator. If the quark is massless, this can happen when 
either p and k are collinear (k \\ p, collinear singularity) or 
when the gluon momentum vanishes (k — > 0, soft singular- 
ity). Thus, a Feynman diagram with all particles massless 
will have soft singularity if it contains an internal gluon 
line attached to two on— shell external quark lines. On the 
other hand, a diagram will contain collinear singularity if 
it has an internal gluon line attached to an on— shell exter- 
nal quark line. It follows then that a diagram containing 
soft singularity contains two collinear singularities at the 
same time, i.e., soft and collinear singularities overlap. 

In view of what has been said above, we conclude that 
the integral lf m is IR finite, and can be calculated in 
D = 4 space— time dimensions, while the rest of integrals 
contain IR divergence (if m and if me collinear divergence, 



I^ mtl , I\ m , and l2 m both collinear and soft divergence), 
and as such have to be evaluated in D = 4 + 2eir {sir > 
0) dimensions. 

For arbitrary D, the integrals |l3|) cannot be ex- 
pressed by elementary functions, but we know that the 
expressions corresponding to these diagrams expanded in 
powers of Sir are of the generic form 

~ 4~ + — + c + 0(e IR ) , 



-IR 



--IR 



with the coefficients A, B, C being complex functions of 
the kinematic invariants. The 1/sjr poles express the IR 
divergence. 

As stated in the Introduction, the IR divergent inte- 
grals If, K e {3m,2m/j, 2me, lm, 0m}, defined through 
Eqs. (0) and (0)— (O) , have been evaluated in Ref. Q 
using the partial differential equation technique. After ex- 
perimenting with various ways to independently derive the 
results of Ref. Q , we have found that the Feynman pa- 
rameter method appears to be the most straightforward 
and satisfactory approach. 






/) 



Fig. 2. One- loop box diagrams with massless internal lines 
but containing 0, 1, 2, 3, and 4 external massive lines; a) the 
four— mass box diagram; b) the three— mass box diagram; c) 
the two— mass box diagram with external masses at adjacent 
corners (legs 3 and 4); d) the two— mass box diagram with 
external masses at diagonally opposite corners (legs 2 and 4); 
e) the one— mass box diagram; f) the massless box diagram. 
Thick lines designate the massive external lines. 

3 Calculation and results 

Using the Feynman parameter method and the method 
of dimensional regularization, in this section we evalu- 
ate the IR divergent scalar one— loop box integrals If , 
K G {3m, 2mh, 2me, lm, 0m}. To accomplish that, we 
start by considering the most complicated of these inte- 
grals, namely, the three— mass box integral i| m . We show 
that, paying due attention to the fact that the limit of 
taking a mass to zero does not necessarily commute with 
the Sir expansion of dimensional regularization, the re- 
sults obtained at the intermediate steps of the calculation 
of the integral /| m can be used to obtain the results for 
the rest of the above integrals. A characteristic feature of 
our calculation is that we keep the causal ie systematically 
through the calculation, so that the results we obtain are 
valid for arbitrary values of the relevant kinematic vari- 
ables. 
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Let us then start with the three— mass box integral 
If 71 . By setting m\ = and D = 4 + 2e IR (e IR > 0) in 
Eq. (^) and eliminating the S— function by performing the 
#4 integration the integral 7| m becomes 



r3m 

A — 



i rp - e IR ) 
(47r) 2 (4:irfi 2 ) eiR J 



dxi 



dx 2 / dx 3 
Ja 



-xix 3 s - x 2 x 3 m 2 . — (1 - x\ — X2 - x 3 ) 



x (x\ m\ + x 2 t + x 3 m 3 ) 



ieir-2 



(14) 



It is a well— known fact that the appropriate choice of 
Feynman parameters is in practice a critical ingredient in 
enabling one to evaluate a complicated Feynman integral 
analytically. There does not appear to be any simple for- 
mula for choosing an optimal set of Feynman parameters 
for a given diagram, as there is generally an enormous set 
of possibilities. 

To proceed with the evaluation of the integral lf m , the 
most suitable set of Feynman parameters turns out to be 
given by O 



x\ = (1 - x) (1 - y), 
X2 =x(l- y), 
x 3 =yz. 



(15) 



The Jacobian corresponding to this transformation of the 
integration variables is y(l — y). Written in terms of the 
new variables, the integral ( J14| ) takes the form 



i r(2- SlR ) 



dxdydz y{l - y)i-y{l - y) 



(47r) 2 (4irfi 2 ) ElR J 
(1 — x)z s + (1 — x) (1 — z) m\ + xz to 2 , + x(l — z) t 

eiR-2 



z (! - z )y 2 m l - ie } 



(16) 



As it is seen from Eq. (|iq), the integration over x is ele- 
mentary and is readily performed. The resulting expres- 
sion is 



r3m 



« [ d dz 1 

2 J V Z z(s-m 2 2 ) + (l-z)(ml-t) 
x| -y(l-y)(z8+(l-z)m: 

- z(l - z)y m 3 - le 

-y(l-y)(zm 2 + {l-z)t) 



z(l 



\ 2 2 
z)y m 3 



SIR — l 



where we have introduced the abbreviation 
i 2r(l-E Jfl ) _ 

(47r) 2 (47T/i 2 ) e -f« 



(17) 



(18) 



Next, by pulling out the factor y SlR 1 from both terms in 
the curly brackets (which is legitimate since y is positive 



in the integration region and the sign of e(> 0) does not 
change), and making use of the following relation: 

(a + b- ie) 6 ^- 1 = f-A- + l] £lR \b - ie) 61 *' 1 
\ b — 16 / 

(a,6eR), (19) 

(which is not self— evident owing to the fact that Ej R is 
not an integer) leads to the integral of the form 



n m = -J dz — -n , 
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z (s — 7713) + (1 — z) (mj — t) 



— z s — (1 — z) to 2 — ie 



eiR — l 



dyy £lR 



sir — l 



i-v 1- 



z(l — z) to 2 + ie 
z s + (1 — z) m| + ie 



■( — z m\ — (1 — z) t — iej 



2 
4 

EIR — 1 



SIR — 1 



SIR— 1 



1-2/1 



z(l — z) TO3 + ie 
zm\ + (1 — z) t + ie 



eiR-l 



(20) 



By noticing that the integral over y stands for the Euler 
integral representation of the hypergeometric function 



iFi(a,b; c;z) 



r(c) 



1 , t b - 1 (i-t) c - b - 1 

dt- K 



r(b)r(c-b)J (\-tzY 

Re c> Re b > ; | arg(l - z) |< tt,(21) 



we obtain the result 

r-l 



i3m 
1 4 



K 



2e IR 



dz 







z s — (1 — z) to 2 — ie 
X2FA 1 — ei R , £z R ; 1 + ei R ; 1 



Z (s — TO 2 .) + (1 — z) (TO4 — i) 

. SIR 1 



z(l — z) to 2 + ie 
zs + (1 — z) m\ + ie 



z m\ — (1 — z) t — ie 



ElR-l 



x 2 Fi 1 — e/B, si R ; 1 + £/_r; 1 



z(l — z) m| + ie 
zto| + (1 — z) t + ie 



(22) 



Next, using the identity (partial fraction decomposition 
1 1 1 



az + bcz-\-d ad — be \a z + b cz + d 



(23) 



and performing a few simple rearrangements, we find that 
the integral under consideration can be written in the fol- 
lowing form: 



i3m 



—(P 3m + Q 3m ) , 

st — m^mX 



(24) 
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where 



P 



3/7/ 



1 



(ml - t) [dz [-z ml-(l-z)t- ie] eiR 
Jo 

z(l 



x 2 Pi( 1 - ei R ,e IR ; 1 + e IR ; 1 



+{m\ - s) [ dz [-z s - (1 - z) m\ - \ef IR 
Jo 



z) m\ + ie 
z m\ + (1 — z) t + ie 

-1 



X2F1I 1 — £/_r, £/_r; 1 + £/_r; 1 



z) m\ + ie 



z s + (1 — z) m 2 + ie 



(25) 



and 



dz 

2 sir J z - z 



(-zm% - (1 - z)t- \ef IR 



x 2 Pi^l - £ir,£ir; 1 + em; 1 
-(-zs - (1 - z)m\ - ie) eiR 

X 2 Pi 1 — Sir, Sir; 1 + Eir\ 1 



z(l — z) m\ + ie 
zm| + (1 — z)t + ie 

z(l — z) m.3 + ie 
z s + (1 — z) m\ + ie 



(26) 



with 



^0 



s + 1 



(27) 

In view of ([m]), it is clear that the integrals if can be 
written in the form 



q = K 9 K {p* 



Q 



K G {3m, 2mh, 2me, lm, Om}, 



with 
ff 3m = 



2rne 



Jlinh 



-,1m 



,0m 



, 995 
SE — 7722 m 4 



(28) 
1 

st 

Not being able to perform the remaining integrations in 
Eqs. (|26| ) and (^H|) analytically, we proceed by expanding 
the integrands in power series in sir and by term by term 
integration. We shall see below that all divergences of the 
integral if are contained in P K , while Q K is completely 
finite. 

It is clear from Eqs. (Il8| ) and (|24|) that to obtain the 
values of the integrals lf~K 6 {3m, 2mh, 2me, lm, Om} 
to order O(e IR ), the evaluation of the integrals P K and 
Q K should be made to the same order. 

It turns out, however, that the number of integrals that 
really need to be evaluted reduces to just four, namely, 
P 3m , given by @, the integrals P 2mh and P 2me , ob- 
tained by setting m| = 0, and m| = in (p5[), respectively, 
and the integral Q 3m given by ( j26| ) . The values of all the 
other integrals appearing in Eqs. ( p8| ) can be derived by 
taking appropriate zero— mass limits. 



3.1 Calculation of the integrals P 3m , p' 2mh , and P 2me 

Let us start by considering the integral P 3m given by Eq. 
( ^5| ) . Making the change z — > 1 — z in the second term on 
the right— hand side in (|25|), one finds that the first and 
second terms are related to each other by the t — ► s and 
m| — > m| interchanges. Therefore, the expression for P 3m 
takes the form 

P 3m = \ [R(t, ml mj) + R(s, m 2 , m 2 3 )] , (30) 

where 

P(a,/3,m|) = 

1 r 1 

dz(f3-a){-zf3-{l-z)a-ie) £lR - 1 



sir Jo 

„ ( z(l — z) m\ + ie 

x 2 Fi 1 - e IR , e IR ; 1 + e 1R \ 1 - a , ,' \ — r-r- 
\ z/5 + (1 — z) a + ie 

(31) 

In order to evaluate this integral, we first apply the linear 
transformation formula for the hypergeometric functions: 



2 Fi(a, b; c; z) 



r(c) r(c-a- b) 
r(c-a)r{c-b) 
x 2 Pi(a, b: 1 + a + b - c; 1 - z) 



r(c) r{a + b- c) 



(l-zf 



—a — b 



r(a)r(b) 

X2-Fi(c — a, c — 6; 1 + c — a — 6; 1 — z) 
c-a-b^±n, I arg(l-z) |<tt. (32) 
With Eq. (p2| ) taken into account, (|3l]) becomes 
i?(a,/3,m§) = 

1 r 1 

= / dz(j3- a)(-z (3-(l-z)a- ie)""" 1 

sir Jo 

r(l + £ IR )r{£m) 



r{2e IR ) 



X 2 Pl 1 — £ir, Sir] 1 — Sir] 



z(l — z) m 2 + ie 
z /3 + (1 - z) a + ie 



—z(l — z) mj — ie 
-z/3 — (1 — z) a — ie 



x 2 Fl 1; 1 + e Jfl ; 



z(l — z) m\ + ie 
z /3 + (1 — z) a + ie 
Assuming, with no loss of generality, that 

I m\ \<\a\, |/3 I a//3>0, 
from which it follows 

z(l — z) + ie 



(33) 
(34) 



z/3 + (1 - z) a + ie 



< 1 



6 [0,1], (35) 



and making use of the series representation of the hyper- 
geometric function 



2 Pi (a, b; C ;z) = ^2 



n=0 



r(a + n)r(b + n)r(c) z n 
r(a)r(b)r(c + n) nf 



z |< 1, 



(36) 
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leads to the result 



R(a, P,ml) 



r(e IR ) 

r(2e IR ) 



(3 + ie 



E 

n=0 



\ [ ' r(i + n) 



x / dz [z(l - z)} 1 



1 - z 1 



(3 + ie 



— (l+n— sir) 



-{-ml-ie) 



tiR r(n + 2e m ) 
r(l + n + ei R ) 



x dz [z(l - z)] 



n+eiR 



/3 + ie 
a + ie 



-(l+n) 



(37) 

Taking into account the integral representation of the hy- 
pergeometric function given by Eq. yields 



R(a,P,m 3 ) = 
(-a-ie)' 



r(e Ia ) 



r(2e IR 

IR r{n + e IR )r(\ + n) 
r{2 + 2n) 



x 2 F 1 [l + n- e m , 1 + n; 2 + 2n; 1 



/3 + ie 
a + ie 

( 2 Y , n r{n + 2e IR )r{l + n + e IR ) 
{ 3 ' r(2 + 2n+2e IR ) 



x 2 Fi 1 + n, 1 + n + sir; 2 + 2n + 2e IR ; 1 



(3 + ie 



a + ie 



(38) 

Although obtained under the restriction expressed by Eq. 
(|34|), the above expression for R(a,(3, m\) can be analyt- 
ically continued for an arbitrary value of (3. The series 
in ( pq ) is convergent because it converges for (3 = (see 
below) and the hypergeometric functions reach the maxi- 
mum at the same point. To retain the convergence of the 
series, we still assume that | ml |< | a \. 

By inspecting the sum on the right— hand side of (|38|), 
we find that, if (3 ^ 0, all terms with the exception of 
the first one (n = 0) are of higher order in Sj R and can 
therefore be omitted. Upon replacing the whole sum by the 
first term, and performing a few simple rearrangements, 
we find 



R(a,(3^0,mi) 



2 , _ r\e IR ) 



1 - 



(3 + ie 



X 2 Fi II -ejij,l;2; 1 
- (-ml - ir^ /R 

x 2 Fi ( 1,1 + em; 



r(2e IR ) \ n + v 

(3 + ie 



(-a - ief IR 



a + ie 



1 



2 (1 + 2e IR ) 



2e 



IR, 



(3 + ie 



a + ie 



+ 0(e IR ). 



(39) 



On the basis of formula ( pT[ ) , one can easily show that the 
first hypergeometric function on the right— hand side of 
the above expression can be written in the form 



F 1 l-e Ifl ,l;2;l 



0- 



(3- 



£l R 



a + ie 

-l 



(40) 



(-a-ie) e 



{-a 



Now, inserting (J40j) into (|39|) leads to the expression 
2r\e IR ) 



R(a,(3,m 3 ) 



r(i 

(-ml - ie) 6lR 



2e IR ) 
sir 



(-a - ie) eia - (-(3 - ie) £,R 
(3 + ie" 



1 - 



2 1 + 2e IR 

x 2 F 1 [ l,l + £ Jfl ;2 + 2£ Jii ;l- 



a + ie 
(3 + ie 



0(e IR ) 



(41) 



which is valid for arbitrary values of a and m§, and for 
(3^0. 

On the other hand, if (3 = 0, the first and second hy- 
pergeometric functions appearing on the right— hand side 
of Eq. ( pl| ) reduce to 

r(2 + 2n)r(e IR ) 



and 



r(l + n)r(l + n + e IR ) 
r(2 + 2n + 2e IR ) r(e IR ) 



r(i 



2e IR )r(\ + n + e IR y 



respectively. This implies that when (3 = 0, all terms in 
the sum are individually divergent and of the same order 
in ei R and, as such, all have to be taken into account, i. e., 
the summation has to be performed explicitly. Therefore, 
as far as the expression for R(a, (3, m§) given by Eq. (|||) 
is concerned, the cases (3 = and (3^0 ought to be 
considered separately. 

By setting (3 = 0, Eq. (|8|) takes the form 



R(a,0,mi) = 

= r 2 (e IR ) 

r(2e IR ) 



(-a - ie) ei 



- (-ml - ier« J2 



1 



1 



mg + ie 



£i R \ a + ie 



n=0 



n + 2s i R 



(42) 



For the purpose of performing the summations in (|42j), we 
note that the series representation of the hypergeometric 
function given by ( |36| ) particularized for a = 1, c = 1 + b 
leads to the formula 



2 F 1 (l,b;l + b;z) = J2 



z\<l, 



(43) 
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which, when taken into account in (|4^), leads to the result 

2r 2 (e IR ) 



R(a,0,mi) 



r(l + 2e IR ) 



(-a-ie) £ ' R 2 F 1 l,e m ;l + e IR 



772.3 + le 
a + ie 



(-m 2 3 - ie) £ ' R 2 F X ( 1, 2e IR ; 1 + 2e IR ; Hll±21 



(44) 

This expression for R(a, 0, m|) is valid for arbitrary values 
of a and 7713 . 

To establish contact with the results obtained in Ref. 
|12| . instead of expanding everything in the expressions 
( [41| ) and ( (4~4| ) to the required order in Ej R , let us, for 
the time being, expand only the occurring hypergeometric 
functions. The relevant expansions are 

z 2 fi (1, 1 + 5; 2 + 25; z) = 



1 + 2S 



l-(l-z) 5 - 2 5 2 U 2 (z) + 0(5 3 ) (45) 



and 



2 F 1 (1,S;1 + S;z) = l-S ln(l - z) - S 2 U 2 (z) + 0{6 3 ) , 

n (46) 

where Li 2 (z) stands for the Euler dilogarithm [|14j defined 
as 

f 1 dt 

Li 2 (z) = - / - In(l-te) . (47) 
Jo 1 

Taking (^5|) into account in (f4l|), we find that 
2r 2 (e IR ) 



R(a,f3,m 3 ) 



r(l + 2e IR ) 



(-a - ie) £lR - (-/3 - ie) £lR 
2 



2 V 3 ; 2 (-a-ie) e « 



£7 fl Li 2 I 1 



+ 



■0( £ 



(48) 



Next, using the expansion (|46|), the expression ( p4[ ) be- 
comes 

R(a,0,mf) = 



-Li 2 1 



TO3 + ie 







\ 7T 2 " 


1 







0(e/«) (49) 



In arriving at (j49|), the relation 



Li 2 (l - z) = -Li 2 (z) - ln(z) ln(l - z) + — , (50) 

6 



has been employed. On the basis of Eqs. (48) and ( |3C| ) we 
find that the integral P 3m is given by 



+ 



(-m 2 - 


ie) £lR - (-ra? - 


- ie) £lR - ( 


1(- 


ml 


_ ie) £ -ffl(-m^ - 


ie) £IR 


2 




(-t - \e) ei * 




1(- 




- ie) eiR (-ml - 


ie) £lR 


2 




(-s - ie) £lR 






Li 2 


( ml+ie\ 
V t + ie J 


+ Li 2 (\ - 



+ 0(e IR ). (51) 
Similarly, combining Eqs. (|48|), ([49|), and (p0|), we obtain 

2m/i _ r2 (£iR) 



P 




(s - ie) £lR + (-t - ie) eji 

(-mj - ie) £lR - (-m 2 , - ie) £lR 
1 (-m| - ie) £lH (-m| - ie) £/R 
2 



(— s - ie) e -f« 
7r 2 . / m l + i £ 

— Ll9 1 ; — 

6 V t + ie 



• Lis 1 



+ 0(e 7fl ). (52) 



We now turn to evaluate the integral p 2me . To this end, we 
set m\ = in Eq. (p5~|). As a result, both hypergeometric 
functions appearing in (B5h reduce to 



2-Fi(l — e/fl, ei R ; 1 + 1) 



r(l + e IR )r(e IR ) 
r(2e IR ) 



, (53) 



making it possible to perform the remaining integration 
analytically. The exact result for the integral P 2me is 



p 2 me 



r 2 (e IR ) 



(s - ie) £IR + (-t - ie) £IR 



r(l + 2e IR ) 
- (-ml - ie) £lR - (-ml - ie) £lR 



(54) 



It is valid for arbitrary values of m\ and m\. Consequently, 
the expressions for the integrals P lm and P 0m are ob- 
tained by setting m\ = and m\ = ml = 0, respectively, 
in®. 

A remark concerning the issue of the zero— mass lim- 
its of the massive integrals P K is in order. By looking 
at the expressions for the integrals P K given above, one 
observes that the limits P 2mh -> P lm , P 3m -> P 2me , 
and P 3m — > P 2mh are not smooth. On the other hand, 
the limits P 2me -> P lm and P lm P 0m are smooth! 
In general, there is no reason for the zero— mass limits 
to be smooth. Namely, the limit of taking a mass to zero 
does not necessarily commute with the l/si R expansion 
of the dimensional regularization, which has been trun- 
cated at (D(e® R ). Note, however, that if we were able to 
evaluate the integral P 3m in (|H) analytically for general 
Ei R , the result thus obtained would suffice to obtain the 
results for the other integrals P 2mft ,F 2me , P lm , and P 0m 
by simply setting m\ = 0, m\ — 0, m\ — m 2 — 0, and 
m\ = m 2 = ml = 0, respectively, and then expanding 
these results to the required order in ei R . 
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3.2 Calculation of the integral Q 3m 

As it is seen from (^), the integral Q 3m is given in terms 
of two hypergeometric functions, both of which can be 
conveniently written in the form 



Fi [ 1 — em, em; 1 + em; 1 



z(l — z) mj + ie 
z a + (1 — z) (3 + ie / 

(a,p)e[(ml,t),(s,ml)}. (55) 

Making use of the transformation formula 

a Fi(o, b; c; z) = (1 - z) c - a - b 2 F 1 (c -a,c-b;c; z), (56) 

the symmetry of 2 Fi with respect to the arguments a and 
b, i.e., 

2 F 1 (a,b;c;z) = 2 F\(b, a; c; z), (57) 

and the integral representation of the hypergeometric func- 
tion (PlD, we can write the hypergeometric function given 
by (|55|) in the form 

z(l — z) m| + ie 



F\ 1 — Em, em; 1 + 1 



—z(l — z) m§ — ie 
-za — (1 — z) [3 — ie 



za+ (1 - z)/3 + ie 



r(i 



-/A 



)r(2 £/fl ) 



x / dyy 28 "- 1 ^-!/)- 6 " 

z(l — z) m\ + ie 



1-y 1- 



z a + (1 — z) [3 + ie 



(58) 



Upon substituting (|58|) into (|26|) , and utilizing the identity 
(E3[) once more, we find that 



r(e m ) 



2r(l-e IR )r(2em) 



dz- 



1 



z - z 



x [-z(l~ z)mi-ie] £ia / dy y 2eiR (1 - y) _ejn 



x { [— «(1 - z) TO3 + z s + (1 - z) ml] 

x \—z s — (1 — z) m\ — ie 

+ y (—z(l - z) m\ + z s + (1 - z) m^)] 1 

- [-2(1 - z) TO3 + z m 2 2 + (1 - z) t] 

x [— z m?, — (1 — z) t — ie 

+ y (-z(l-z)m2 + zm2 + ( 1 _ z ) < )]- 1 | 



(59) 



Since we are interested in obtaining the value of Q 3m to 
0(s° m ), the fact that the expansion of the prefactor in 
the above expression is of the form 1 + O(em), allows us 
to set Em = in the integrand in (p9|). As a result, the 
expression for Q 3m reduces to 



Q 3 ™= fdz^— f 1 
Jo z ~ z o Jo 



dy 



(60) 



{ [-z(l - z) ml + z s + (1 - z) : 



x [— z s — (1 — z) 771,4 — ie 
+ y (—z(l — z) m\ + zs + (1 — z) mf)] 
- [-z(l- z)ml + zml + (l-z)t] 
x [— z ml — (1 — z) f — ie 

+ y (-z(l -z)m^ + z ml + (1 - z) t)] _1 }+ 0(e/ fl ). 
Performing the y integration, we find 



Q 3m = 



■In 



z(t — m|) — i — ie 



o z - z \z{m\-s)-m 



ie 



0(e/fl) 



( 61 ) 

To carry out the remaining integration, it is important 
to note that the residue of the integrand at the pole zo 
is zero, and the logarithm does not cross the cut. This 
fact allows us to make a few simple transformations of the 
integrand. 

Thus, upon the substitution z — » z + Zq, the decompo- 
sition 

pi r-zo A ~ z Q 

dz —* — / dz + dz , 

10 Jo Jo 

followed by a change of the variable z — > — z zo in the first, 
and z — > z(l — zo) in the second term, the integral Q 3m 
can be written down as 



0/ 3m ) 

c3m\~\ 



where we have introduced the abbreviation 





\z 








{Mi 


Jo 


z 


+ ln 


[1 


-z(l 


-In 


[1 


-z(l 


-In 


[1 


-z(l 



/•3m 



s + t — m. 



2™2 



St — TO2 m 4 



(63) 



Expressed in terms of the Euler dilogarithm, the final re- 
sult for the integral Q 3m takes the form 



Q 3m = Li 2 1 - (s + ie)/ 



3 m 



-Li 2 



1 - (t + ie)/ 



3 m 



(64) 



-Li? 



l-(mi+ie)/ 



3rfi 



-Li? 



l-(mi + ie)/ Jm +0( £jii ) 



The expression for Q 3m does not depend on m|, and is 
valid for arbitrary values of m\ and m\. A consequence 
of this is that the expression for the integral Q 2mh to 
the same order in em can simply be obtained by setting 
m\ = in (U). Note, however, that, strictly speaking, the 
expressions for Q 2me , Q lm , and Q 0m cannot be obtained 
by taking appropriate zero— mass limits of the same ex- 
pression. Namely, the expression (^9|) has been derived 
assuming that m\ ^ 0. Therefore, the Em expansion of 
([59]) is not justified in the m§ — > limit. In order to obtain 
the integrals Q 2me ,Q} m , and Q 0m , we proceed as follows. 
We return to Eq. (|26| ) and set m\ = 0. As a result, Eq. 
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!) 



I) reduces to 



)2me 



r\e IR ) 



1 



+ 2 Li 2 [ I 

[(-zm 2 -(l- z)t-ie) eiR +2Li 2 



m\ + ie 



2r(2ei R )J z-z 
- (-zs -(l-z)m\- ie) eiH ]. 

Expanding this expression into power series in si R 



t + ie 
l-( s + ie)/ 3m 



2Li 2 1 



mf + ie 



2 Li, 



1 - (t + ie) / 



3 m 



(65) 



-2Li 2 
-2Li 2 



Q2me = 



' dz , ( z(t — m 2 ) — i — ie 



l-(m 2 2 +ie)f m 
l-(m 2 +ie)f m V l O! < ;; ,i 



(69) 



■In 



o z- z 



z(m\ 



0(e m ), 
(66) 

and comparing it with the expansion of Q 3m given by 
Eq. (61), we find that the expansions for Q 3m and Q 2me 
coincide to the required order in ei R . This being the case, 
all other intgerals Q K can be derived from the integral 
Q3m gj ven by gq ( j54| ) by taking appropriate zero— mass 
limits. 



3.3 Results 

Having obtained, in the preceding subsections, the closed 
form expressions for the integrals P K (s, t; {rrvf }) and 
Q K (s, t; {m| }) to order O(e IR ), we are now in a posi- 
tion to write down explicit expressions for the integrals 

Jf(M;{m?}). 

Before proceeding, it is convenient to introduce the 
functions 



the adjacent ("hard") two— mass scalar box integral is 

sir 



j2mh 



(s,t;m 3 ,m 4 ) 



-t-ie 
4nfi 2 



l rp 



(4tt) 2 st 



— ms — ie 



'IR 



Ann 2 



—mj — ie 



+ 



'IR 



Ann 2 



-mi — ie 



+ 2Li 2 I i - 
+ 2Li 2 

-2Li 2 



mj + ie 



-2Li 2 1- 



m 3 + ie 



1 - (s + ie) f 2mh 
1 - (ml + ie) f 2mh 



2Li 2 

\ + 0(e m ) 



t + ie 
1 - (t + ie) f 2mh 



(70) 



the opposite ("easy") two— mass scalar box integral is 



r n = f 

f 2mh = f 



2me 



s + t — m 



Ii me ( Sl t;mlmj) 



rr 



st — m\m\ 



(Air) 2 st — m\ m 2 



s + t- m\ 

si 
s + 1 

~sT ' 



'IR 



4:TTfi 2 



-t-ie 

4:TTfi 2 



4 

SIR 



— ms — ie 



4ir[i 2 



and the notation 



rr 



r(l-e IR )r 2 (l + e IR ) 
r(l + 2e IR ) 



(67) 



(68) 



-mj — ie 



4lT{i 2 



2 Li 2 



1 - (s + ie) f 



2 m c 



+ 2Li 2 
-2Li 2 



1 - (t + ie) f 2me 
l-(m 2 + ie)/ 2me 



- 2 Li 9 



l-(ml + ie)f 2me 
0(e IR ) . (71) 



for the r function prefactor. 

Now, on the basis of Eqs. (|28|) , (||), 
(|54|), we find that 

the three— mass scalar box integral is 



/| m (s,t; m 2 ,mimi) 



r r 



-IR 



-s — ie 

47T^ 2 

Sir 



(Air) 2 st — m 2 m\ 



-t - ie 
Air^i 2 



")T|) (F2h and ^ n essen tial feature of the above expression for 7| m e is 
' that it is well behaved in the m\ — > 0, m\ — > limits. A 

consequence of this is that it contains both the one— mass 
scalar box integral and the massless scalar box integral. 
Thus, setting m\ = in (|7|) and making use of the fact 
that Li 2 (l)=7r 2 /6, we find that the one— mass box scalar 
integral is 



47T/1 2 



47T/i 2 



+ 



'IR. 



2 ■ \ £ IR 

-m% — ie 



47T^ 2 

,2 



1 / 2 ' \ sir / 2 • \ sir , 

1 / — m| — ie \ / — m A — ie \ ( —s — ie 



'IR 



4-KfJ, 2 



4lT{i 2 



4ir[i 2 



-t - ie 



Il m {s,t;mi) 



1 rp 



(4tt) 2 st 



'IR 



-s — ie 
47r,u 2 



+ 



-t - ie 



4lT{i 2 



2 ■ \ Sir 

-m, — ie 



47T^ 2 



47T/1 2 



47r^ 2 



4-np 2 



2Li 2 
2Li 2 



1 - (s + ie) f 1 



2 Li 



1 - (ml + ie) f lm 



1 - (t + ie) f lm 
\+0(e IR ). (72) 
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Finally, setting m\ — in ([72]), we find that the massless 
scalar box integral is given by 



4 (S ' i} -(4^ st\ej R 
2Li 5 



4:TTfl 2 



2 Li, 



1 - (t + ie) f 



On) 



1 - (a + ie) f m 



2y +0(e/K). (73) 



The above expressions for the one— loop IR divergent scalar 
box integrals I^(s,t; {mf}) constitute the main result of 
this paper. It is important to emphasize that, owing to 
the fact that we have kept the "causal" ie systematically 
throughout the calculation, these expressions are valid for 
arbitrary values of the relevant kinematic variables: exter- 
nal masses mf (i = 2, 3, 4) and the Mandelstam variables 
s and t. 

As stated in the Introduction, the integrals if {s, t, {mf }) 
have been evaluated in Ref. with the help of the partial 
differential equation technique. The calculation has been 
performed in the Euclidean region, where all kinematic 
variables are negative, i.e., 



s, t < 0, m 



2- 



'3' 



< 



(74) 



and the results thus obtained have been analytically con- 
tinued to the positive values of the kinematic variables 
(physical region) by applying the following replacements: 



t — » t + ie , m. 



(75) 



In order to facilitate the comparison of our results with 
those of Ref. |T^ ] , we have written our results in the same 
form in which they have been presented in Ref. fl2]| . A 



glance at the expressions (69-|73J) reveals that they all 



have the same general form, namely, 



l£(s,t;{mj}) = 



(47T) 



r r g 



G K (s,t;s IR ;{mf}) 



~-IR 



arrived at the following conclusions: First, the results are 
in complete agreement in the Euclidean region. Second, for 
the integrals Ij mh (s, t; m|, m|), l| m (s, t; mf), and 7f"(s, t), 
we have found agreement for arbitrary values of the kine- 
matic variables. Third, the results for the integrals 
Il m (s,t; ,m\,m\,mf) and I% me (s,t;m%,mf) do not agree 
outside the Euclidean region. 

The reason for this disagreement is that the analyti- 
cal continuation from the Euclidean to the physical region 
as given by Eq. ([7q) is not well defined for all terms ap- 



+ H K (s,t;{mj})} + 0(e IR ) 
K G {3m, 2mh, 2me, lm, 0m}. (76) 

The IR divergences (both soft and collincar) of the inte- 
grals are contained in the first term within the square 
brackets, while the second term is finite. The function 
G K (s, t; Sir; {mf }) is represented by a sum of powerlike 
terms, it depends on Ejr and is finite in the Ejr — > limit. 
As for the function H K {s,t; {mf}), it is given in terms of 
dilogarithmic functions and constants. 

Comparing our results with the corresponding ones of 
Ref. p2| , we find that the expressions for G K (s, t; cir; {mf}) 
are in agreement. On the other hand, the correspond- 
ing expressions for the terms H K (s, t; {mf}) are of dif- 
ferent form. Proving the equivalence of our results for 
the integrals I± (s, t; {mf }) with those of Ref. Jljj then 
amounts to showing that the expressions for the terms 
H K (s, t; {mf}) agree numerically. By doing this, we have 



pearing in the expressions for the integrals I^(s,t;mf^ 



point 

applying the replacements (|75|), no cut will be hit by the 
powerlike terms, logarithms, and the dilogarithms with 
a single ratio of the kinematical variables. In addition 
to this kind of terms, the expressions for the integrals 
ll m {s, t; , ml, ml, ml) and /| me (s, t; m\, mf), given in Ref. 
[[L2| , contain terms of the form 



Li 2 1 



st 



(77) 



i.e., the dilogarithms of a product of ratios of the kine- 
matic variables. This type of term requires special care. 
In order to avoid crossing a cut, in this case one has to 
make the following replacements: 



Li 2 1 



Li 2 1 



ie m.4 



st 



-V[ 2 , . , ,\. Unl 1 

s + ie t + ie 



where the function r/ is defined as 



t + ie 



m\ + ie m\ + ie 
s + ie t + ie 



(78) 



T)(x, y) = \n(xy) - \n(x) - \n(y), 



(79) 



and arises from the possibility that x, y, and xy are not 
on the same Riemann sheet. Choosing the principal value 
of the logarithm such that the cut lies along the negative 
real axis, Eq. (|79|) can be written as follows: 



(80) 



rj{x,y) — 27ri{#(— Im x)0{— Im y)(9(Im xy) 
— 0(lm x)6(lvn y)0(— Im xy)} . 



If the terms of the form given 
continued in accordance with 



and if me 



If 

t;m%,; 



( [77|) are analytically 
we find that the re- 



sults of Ref. ||1 2| for the integrals Il m {s,t;m1, 



2 2 \ 

m%,mi) 



(s, t; m^, mf) are numerically equivalent to the 
corresponding results obtained in this paper for arbitrary 
values of kinematic variables. 

Having thus numerically established the equivalence of 
the two sets of results for the integrals I^(s,t; {mf}), our 
next task is to demonstrate this equivalence analytically, 
i.e., explicitly. There are, in principle, two approaches to 
accomplish this. The first approach consists in applying a 
series of the Hill identities |l4|] (relating the dilogarithms of 
different arguments]to the final expression for the integral 
Q3m gj ven by jrjq ( |(34| ) t with the aim to express it in terms 
of the dilogarithms occuring in the final expression for the 
integral lf m {s, t; {mf}) in Ref. However, because of 
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the presence of ie in the arguments of the dilogarithms 
in Eq. (Q), this turns out to be extremely messy. In the 
second approach, one tries to achieve the same by recalcu- 
lating the integral in (|6l| ) with an appropriate change of 
the integration variable. After a lengthy calculation, de- 
tails of which are presented in Appendix A, we have been 
able to show that 



)3m 



= -Li 2 1 - 



- Li 2 
+ Li 2 
+ ry 



nr, 



Li 2 1 



lr 



s + ie 



Li 2 1 



nr. 



t + ie 
ic" 



ie m\ 



t + ie 



t + ie 
2 I s + ie 



m\ + ie m 2 + ie 



s + ie 



t 



In 1 



t + ie 
ie m% - 



t + ie 



+0(e IR ) 



(81) 



Although different in form, one can readily prove that 
this expression for Q 3m is numerically equivalent to that 
given in Eq. ([34]) for arbitrary values of kinematic vari- 
ables. All other integrals Q K ,K G {2mh,2me,lm,0m} 
can be derived by taking appropriate zero— mass limits. 
Combining the expressions thus obtained for Q K with ex- 
pressions ([3l]), (p2[), and ( |54| ) for P K , we arrive at the 
final expressions for the integrals l¥ (a, t; {raj }) which are 
in agreement with those of Ref. [|12| if provided the latter 
are correctly analytically continued outside the Euclidean 
region. 



4 Conclusion 

Using the Feynman parameter method, we have calcu- 
lated in an elegant manner a set of one— loop box scalar 
integrals with massless internal lines, but containing 0, 1, 
2, or 3 nonzero external masses. To treat IR divergences 
(both soft and collinear), the dimensional regularization 
method has been employed. We have kept the causal ie 
systematically throughout the calculation. Consequently, 
the results for these integrals, which appear in the pro- 
cess of evaluating one— loop (N > 5)— point integrals and 
in subdiagrams in QCD loop calculations, have been ob- 
tained for arbitrary values of the kinematic variables and 
represent the extension of the results of Ref. outside 
the Euclidean region. 

This work was supported by the Ministry of Science and Tech- 
nology of the Republic of Croatia under Contract No. 00980102. 



analytical continuation for the terms of the form given 
in ([77|), which appear in the expressions for the integrals 
^"(s^jm^mim^) and /f me (s, t; m^, 777,4) of Ref. JlJ, is 
given by Eq. (|78|). 

To accomplish that, we return to the integral Q 3m 
given in Eq. (^Tj), in which, for convenience, we replace 
the integration variable z by y. Passing to the new inte- 
gration variable given by 



1 



y(t -mp-t-ie 
y{m\ - s) — mf - ie 



the integral Q 3m becomes a line integral 

rZ2 dz s + 1 — m\ — to| 

Zl z z{m\ — s) + s + t — m\ - 
+ 0(e IR ) , 



Q 



3 m 



(82) 



ln(l - z) 



(83) 



with the path of integration followed from z\ to z 2 , where 



Zl 



1 



t + ie 



Z2 



1 



my, 



ie 



s + ie 



(84) 



Applying the partial fraction decomposition (|23j), this in- 
tegral can be represented as 



where 



and 



Q 3m = I + J +0(e IR ) 
" Z2 dz 



ln(l - z) 



(85) 
(86) 



Zl 



J 



Z'2 



dz 



z{m 2 z — s) + s + t 



m% — m% 



Tn(l - z). 



(87) 

Let us now consider these two integrals in turn. 

The integrand in ( |S6| ) has a first— order pole at z = 
and the logarithmic branch cut extending from 1 to 00. 
It follows from Eq. (M) that, depending on the values of 
the parameters s, t, mg, ml, it might happen that the line 
connecting the points z\ and z 2 crosses the real axis. In 
this case, as it can be seen from Eq. (|82]), the crossing 
occurs at the point z = — the pole of the integrand. 
Observe, however, that the residue of the integrand at 
z = is zero. A consequence of this is that, regardless of 
whether the line connecting z\ and z 2 crosses the real axis 
or not, we are allowed to assume that it passes through 
the point z = 0. Consequently, the line between zi and z 2 
can be decomposed into two segments: one connecting z\ 
and 0, and the other connecting and z 2 . In view of this, 
the integral / can be rewritten in the form 



/ = 



* 2 de 



ln(l - z) 



z i dz 



ln(l - z) 



Appendix A 

In this Appendix we analytically demonstrate that for Eu- 
clidean kinematics our results for the intgerals I± (s, t; {mf }) 
are in a gree ment with the corresponding results obtained 
in Ref. p2| . As a byproduct, we prove that the correct 



By changing the integration variables z — > z z 2 and z — > 
z zi in the first and second integral, respectively, and tak- 
ing formula (f47|) into account, we obtain the result 



I = -Li 2 1 - 



TO2 + ie 
s + ie 



Li *"-sn!y- (89) 
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Applying the transformation 

Li 2 (l - ij = -Li 3 (l ~y)- \ In 2 V (90) 



to the second term on the right— hand side in (B9J), the 
final expression for the integral I is found to be 



/ = -Li 2 1 - 



to 2 + ie 



s + ie 

.2 

- ^ In , 

2 V i + ie 



- Li 2 1 - 



t + ie 



(91) 



Turning now to the integral J in (jS7|), we switch back 
to the old integration variable as a result of which the 
integral takes the form 



J = J i + J 2 + J 3 



(92) 



where 



Jl = ~ / d ^r7-2 — is — —t 



y {m\ — s) — m\ — ie 
In [y(m\ — s) — m\- ie] , 



(93) 

= I dy __,__ 2 m \~ a _ 2 . In (-* - ie) , (94) 



o y (ml - s) — m\ — ie 



and 



o y(m\ - s)-m\- ie 



t — m? 



The integrals Ji and J 2 are elementary, and are readily 
evaluated. The results are 

Jl = ~ [ ln 2 (-s - ie) - ln 2 (-m2 _ ie ) ] , ( 96 ) 



J 2 = In (— t — ie) In [ — 5— 

\m| + ie 



(97) 



In order to evaluate the integral J 3 , we proceed by adding 
and subtracting the following expression: 

:i V \ t + ie J J y 

dy 1 , ln(l - y) . 

m| + ie 

2/ 2 

mj — s 

After a simple algebraic reduction, the integral J3 can be 
represented in the form 



^3 — ^3,1 + ^3,2 



(98) 



where 



v \ t + le 



d y -2 - ■ H 1 - v) 

TO| + ie 



v 



y -^2 — 

mi — s 



= -Li 2 1 - 



ms + ie \ 7r 



, , --Li 2 (l- W4 + i£ 
t + ie J 6 \ s + ie 



(99) 



and 



^3.2 — 



dy 



In 



1 



V 



t — mt 

y , , ■ - / 

i + ie / 
(. 10 °) 

Next, consider the integral J 3 . 2 . Introducing a new inte- 
gration variable 



i-y 



t 



(101) 



t + ie 



and the notation 

_ (m\ + ie) (777.4 + ie ) 
(s + ie)(< + ie) 

the integral can be cast into the form 



■Ha — 



1 dz 



ln(l -z) + 



dz- 



1 



(102) 



In z. (103) 



1 - a 



Notice that, in a ccor dance with (101), the integral J 3i2 
is given by Eq. (103) as a line integral in the complex 
z— plane with the integration path having only the end— po- 
ints on the real axis at z — and 2=1. Adding and 
subtracting the integral of the form 



dz- 



1 



1 - a 



ln(l - a) 



in the second integral in (|103j) allows us to write the inte- 
gral J 3j2 as 



•h,2 = [— ln(l - z) + [ 
Jo z Jo 



dz- 



1 



-[In 2 + ln(l - a)] 



f 

Jo 



dz- 



ln(l - a) 



(104) 



1 



It should be observed that the residue of the first integral 
on the right— hand side in (101) at the pole z = is equal 
to zero. The same is true for the second integral at the pole 
z = 1/(1 — a). Therefore, the integration path in both of 
these integrals can be taken to follow the real axis from 
to 1. After evaluating the first two integrals in ( 104 ) 
and passing to the old integration variable in the third 
integral, we arrive at the following expression for J 3j2 : 



■h,2 = -— + Li 2 (l - a) +ln(a)ln(l - o) + ln(l - a) 
6 



< I dy 

Jo [y\ m 



st — m\m\ + ie(s + t — m 2 — m|) 



'AW 



t - ie 



(105) 
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Carrying out the remaining integration, we obtain the fi- 
nal expression 



J; 



3.2 



6 



-Li 2 (l— a)+r) 



i\ + ie m% + ie 



t + ie 



ln(l-a) , 
(106) 



where the function r)( x, y) is defined by Eq. (|79|). Now, 
substituting ( |99| ) and (10C) into (p^), we find the integral 
J3 to be given by 



■h = -Li2 
+ Li 2 
+ r] 



1 - 



1 



- Li 2 1 - 



mi 2 + ie 
t + ie 



?n| + ie m\ + ie 



s + ie t + ie 
m\ + ie m 2 , + ie 



s + ie t + ie 



In 1 



m/ + ie mi + ie 



ie t + ie 



(107) 



Finally, having obtained all the necessary ingredients, we 
now combine them to obtain (81), which is the desired 
result. 



Appendix B 

By expanding the power like terms appearing in Eqs. (p9|) — 
(f73|), and the prefactor rp defined by d68|), we find that 
the integrals under consideration can be written in the 
generic form 



J 4 -J4^ 9 



K 



A 



K 



+ 



D 



K 



-IR 



£lR 



K G {3m, 2mh, 2me, lm, 0m} 



(108) 



The functions g K appearing above are defined by (p9|). For 
convenience, the finite parts have been decomposed into 
two terms where the term originates from the expan- 
sion of the product of the rp prefactor with the powerlikc 
terms. 

The double— pole parts A K , the single— pole parts B K , 
as well as the finite parts and of the individual 
integrals are listed below. 



3m 



A 



, 
in 



s + ie 
mi + ie 



In 



Of 



3m 



+ 



1 



1 

+ 2 



7e + In 
7-E + In 
IE + In 
7E + In 
In 



Ann 2 



t + ie 
m\ + ie 



7e 



In 



-t - ie 

AlTfJ, 2 



4lT{i 2 



IE 



In 



AlTfi 2 



Anii 2 



47T/1 2 



IE 



—nit — ie 



In 
In 



to| + ie 
s + ie 



t + ie 



-1 2 



£i3m 



2Li 2 I I 
2Li 2 f 



t + ie 
l-(.s + ie)/ 3m 



2 Li, 1 



- 2Li 2 

B 2mh = 1e + m 



1 - (mj + ie) f 



2 Li, 



-2Li 2 



m\ + ie 
s + ie 

3 m 



I- (t + ie) f 

l-{m\ + ie)f m . 

(109) 



-t - ie 



Ann 2 



In 



In 



t + ie 
mi + ie 



£t2mh 



7e + In 



-s — ie 

47T/1 2 



IE + In 



—m% — ie 



4lT{i 2 



JE + In 



Airn 2 



m 2 A + ie 

7E + In 
7e + In 
In 



/ ml 


*■*)] 


V «4 


le J. 



-t - ie 

47T/1 2 

—ml — ie 

AlTfl 2 

2 



7T 
12 



G 



2mh 



2 Li 2 1 



2Li 2 



,s + ie 

2mh 



l-(s + ie) f 



- 2 Li 2 I I 
2Li 2 



2Li 2 l_( m 2+i c )/ 



2mfi. 



< + ie 
1 - (t + ie) f 2mh 

(110) 



A 2me = Q ^ 
B 2me = 2 hi 



m 2 + ie 



2 In 



2me 



7e + m 



47r/i 2 



7-E + In 



C 2 2me = 2 Li 2 



— tos — ie 



2 



47T^ 2 

1 - (s + ie) f me 



t + ie 
m\ + ie 

IE + In 

7-E + In 
f 2Li 2 



-t - ie 

4tT[I 2 



— mi — ie 



-2 Li, 



1 - (ro| + ie) f 



-2Li 2 



47T/K 2 

1 - (t + ie) f 2me 

l-( m \ + ie)f™ . 

(Ill) 



A Lm 



2 7£; + 2 In 



47T^ 2 



2 In 



- ie 

47T/7 2 



2 In 



7e + In 



7E + In 



-mi — ie 

47T^ 2 

— s — ie 



47T^ 2 



47T^ 2 



7e + In 



- ie 

AlTfl 2 



2 2 



G 



2Li 2 l-(s + ie)/ lm +2Li 2 l-(i + ie)/ 



14 



G. Duplancic, B. Nizic: Dimensionally regulated one— loop box scalar integrals with massless internal lines 



aQtu 



Win 



- 2Li 2 

4 , 
47b 

1e - 
2Li 2 

2* 
3 



1 - (m 2 + ie) f 1 



2 In 



7T 

y 



(112) 



4:TTfl 2 

s — ie 



2 In 



-t - ie 



47T^ 2 

( S + ie)/ 0m 



7b + In 
f 2Li 2 



4:TTfl 2 
-t 



TV 

y 



i 



(t + ie) f m 

(113) 

0.5722 is the Euler con- 



In the above expressions, 
stant. 

The pole and finite parts for the integrals I± (s, t; m?) 
obtained in Ref . Jl2| , but with the correct analytical con- 
tinuation outside the Euclidean region for the integrals 
/f m (s, i; m|, m§, m|) and 7f me (s, t; m§, m 2 ), can be ob- 
tained simply by replacing the terms (s,t;m 2 ) given 
above by the corresponding values listed below: 



In 2 



t + ie 

2 i 



2Li 2 
2Li 2 
2 V 



C 2mh = _ ln 2 



- 2 Li 2 1 - 



2 Li 2 1 



t + ie 



ie 



s + ie 
ie m| - 



< + ie 
ie' 



t + ie 



s + ie 
t + ie 
m 2 + i< 
t + ie 



In 1 - 



s + ie t + ie 



(114) 



- 2 Li 2 1 - 



t + ie 



(115) 



Cl me = - In 2 



2Li 2 
2 Li 2 



t + 




m| 


+ ie 


s - 


hie 


ml 


+ ie 



+ 2 Li 2 1 



2r, 



ie 

f ie m 



2 Li 2 
2 Li 2 



t + ie 



t + ie 



ie 



ie 



i + ie 

ie' 



t + ie 



In 1 - 



m 2 + ie ttt.4 + ie 
s + ie t + ie 



(116) 



2 Li, 1 



in > 



In 2 



s + ie 

s + ie 
t + ie 



2 Lia 1 



t + ie 



(H7) 
(118) 



As it is readily seen from the above expressions, the in- 
tegrals 1^ (s,t; {m 2 }) are real in the Euclidean region . 
Outside of this region, however, the integrals acquire an 
imaginary part. Being given in terms of logarithms and 
dilogarithms, their imaginary parts can be easily deter- 
mined. 

Thus, with the usual definition of the logarithms on 
the branch cut — oo < z < 0, one has 



ln(y ± ie) = In \y\ ± m0(-y) 



(119) 



Next, as it is seen from Eq. ((47|), the function Li 2 (y) de- 
velops an imaginary part for y > 1, and 



1 



W,U-An±y ) = -U 1 [^ -- m 2 y+^-, (120) 



Im Li 2 (y ± ie) = ±7rln y. 



(121) 
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